Magic radioactivity of 252 Cf 
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We show that the sharp maximum corresponding to 107 Mo in the fragment distribution of the 
252 Cf cold fission is actually a Sn-like radioactivity, similar to other decay processes in which magic 
nuclei are involved, namely a-decay and heavy cluster emission, also called Pb-like radioactivity. It 
turns out that the mass asymmetry degree of freedom has a key role in connecting initial Sn with 
the final Mo isotopes along the fission path. We suppose the cold rearrangement of nucleons within 
the framework of the two center shell model, in order to compute the cold valleys in the charge 
equilibrated fragmentation potential. The fission yields are estimated by using the semiclassical 
penetration approach. We consider five degrees of freedom, namely the inter-fragment distance, the 
shapes of fragments, the neck parameter and mass asymmetry. We found an isomeric minimum 
between the internal and external barriers. It turns out that the inner cold valley of the total 
potential energy is connected to the double magic isotope 132 Sn. 



S 1 



(N 
> 

t> 
in 

co 
o 
o 



The spontaneous cold rearrangement of nucleons from 
an initial nucleus to two final fragments corresponds to 
the most favorable path of the cold splitting in the po- 
tential energy surface This path, called cold valley, 
is related to the magicity of one or both fission frag- 
ments. Thus, a-decay is connected to the cold valley of 
the double magic nucleus 4 He. The cold valley of the 
double magic 208 Pb is responsible for various heavy clus- 
ter decays, in which C, 0, Ne, Mg and Si are emitted @ 
and this is the reason why they are also called "magic- 
radioactivities" . 

On the other hand the production of superheavy el- 
ements is connected with the inverse fusion process Q, 
involving double magic nuclei 208 Pb and 48 Ca 0, ||. 

Two main synthesis procedures of superheavy nuclei 
based on fusion reactions have been experimentally used: 

(a) the cold-fusion at GSI Darmstadt with either 208 Pb 
or 209 Bi target leading to a small excitation energy of the 
compound nucleus followed by a single neutron evapora- 
tion [1,0], and 

(b) the hot fusion with 48 Ca projectile at JINR Dubna, 
in which the compound nucleus is very excited hence 
more neutrons are evaporated [!, H[ . 

Almost two decades ago, systematic measurements 
we re p erformed to determine cold fission yields of 252 Cf 
[gLfTof. The aim of this work is to analyze these data and 
to show that the cold fission of 252 Cf is strongly connected 
with the cold valley of the double magic isotope 132 Sn, 
although the experimental cold yields have a maximum 
corresponding to a different charge number. 

In the past the cold fission of 252 Cf was inve stig ated 
within the double folding potential method fill . Il2l | em- 
phasizing the role of higher deformations. In Ref fliij . 
based on a macroscopic model by determining the tip 
distances for the exit point from the barrier for ground 
state deformed fragments, it was predicted that the ma- 
jor contribution in the yield distribution corresponds to 



the light fission fragment A2 «100, contradicting exper- 
imental data showing a peak at A2 = 107. In this study 
we explain the root of this discrepancy. 

We extend the analysis performed in Ref. 14] to a 
more reliable microscopic approach to estimate the fission 
barrier, based on a new version of the Super Asymmetric 
Two Center Shell Model 15]. This approach was already 
used to describe the fusion/fission of some superheavy 
elements [lj| and of the dynamical effects in fission [l5|, 

In this work the nuclear shape parametrization is ob- 
tained by smoothly joining two spheroids with a third 
surface, given by the rotation of a circle around the sym- 
metry axis. This parametrization is characterized by 
5 degrees of freedom, namely the mass asymmetry, the 
elongation R — z 2 — Z\ given by the distance between the 
centers of the nascent fragments, the two deformations 
of the nascent fr agments associated to the eccentricities 
£ i = vb^iK (* = 1>2), and the necking characterized 
by the curvature C — S/ R3 of the median surface and 
a asymmetry term given by the ratio of the semi-axis 
r\ = ai/a2. In this parametrization the values 5=1 and 
S=-l stand for necked and swollen shapes in the me- 
dian surface, respectively. Within these generalized co- 
ordinates, for necked shapes it is possible to estimate the 
mass of the nascent fragment A2 during the deformation. 

The penetrability, corresponding to some binary par- 
tition, defines the isotopic yield. This quantity can be 
estimated by using the semiclassical integral [l8[ 
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between internal and external turning points along a fis- 
sion path. Two ingredients are mandatory in order to 
evaluate the action integral: the inertial parameter M 
and the fragmentation deformation energy V (we will call 
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FIG. 1: Deformation energy V computed within the 
microscopic-macroscopic method with respect C (necking co- 
ordinate) and R (elongation). In the lower plot, the mini- 
mal action trajectory is also plotted. The step between two 
equipotential lines is 1 MeV. 



it simply deformation energy). In Eq. ([T]), the inertia M 
is considered along a given fission path where the main 
coordinate is the elongation R. Therefore, for this trajec- 
tory in the configuration space, the dependencies of all 
generalized coordinates versus the main coordinate R are 
known and the quantities V and M depend only on the 
variables A-i and R. The problem of finding the fission 
trajectory will be treated below. For a fixed combination 
A = A\ + A2 the deformation energy has a minimum 




at the charge equilibration point Z2, which we will not 
mention in the following. By definition it is defined as 
follows 

FIG. 3: Deformation energy minimized statically with respect 
V(A2, R) = Vn(A2, R) + Vc(A2, R) — Q (2) the eccentricity of the second fragment £2 and the mass asym- 

metry parameter r\ as function of R and A2 in the second 

where Vn(A2,R) is the nuclear and Vc(A2,R) Coulomb saddle region. 

inter-fragment potential. We also introduced the Q-value 

in terms of the difference between binding energies of the 
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FIG. 4: Deformation energy V computed within the 
microscopic-macroscopic model for different binary partitions 
with respect A2 and the elongation R. In the lower plot, the 
inertia are plotted. 



parent and the sum of emitted fragments, i.e. 

Q = B(Z 1 , A-l) + B(Z 2 , A 2 ) - B(Z,A) . 



(3) 



For deformed nuclei, due to the fact that the largest emis- 
sion probability corresponds to the lowest barrier, the 
deformation potential decreases in the direction of the 
largest fragment radius. 



3000 



2500 



Z5 

6 
>- 



2000 



1500 



1000 



500 



iiv 




116 




i> ' 

/ f 

— / 1 




/ 1 

/ 1 
/ ' 

4>0© 


1 
1 

©3 

Q 

© 




0-0 v 
© b 


95 100 105 


110 115 120 

A 2 



FIG. 5: Experimental yields in arbitrary units (dashed line), 
compared with renormalized theoretical penetrabilities calcu- 
lated within the microscopic-macroscopic model (solid line) 
with respect to A2. 



The deformation energy of the nuclear system is the 
sum between the liquid drop energy Vr.p and the shell 
effects 6E, including pairing corrections i.e. 



V = V LD + SE- V 



(4) 



The energy of the parent nucleus Vq is used as a refer- 
ence value, so that in the ground state configuration the 
deformation energy is zero and asymptotically, for two 
separated fragments at infinity, the deformation energy 
reaches the minus sum of energies of emitted fragments. 
Thus, Eq. ((H) coincides with the definition Eq. fl2j. 

The macroscopic energy is obtained within the frame- 
work of the Yukawa-plus-exponential model [20| . ex- 
tended for binary systems with different charge densities 
[2]], [22)]. The Strutinsky prescriptions [23| were computed 
on the basis of a new version of the superasymmetric two- 
center shell model. This version solves a Woods-Saxon 
potential in terms of the two-center prescriptions as de- 
tailed in Ref. 15]. The inertial parameter M is com- 
puted in the framework of the cranking model 0, [25| . 
We considered only cold fission process. Consequently 
the deformations of the final nuclei are given by their 
ground state values of Ref. '26i\ . 

For comparison with experimental data, the maximal 
values of the independent yields for a maximum excita- 
tion energy of 7 MeV were selected from Ref. @. The 
selected channels address binary partitions characterized 
by the following light fragments: 95 Rb, 96 Rb, 97 Sr, 98 Sr, 



3: 95 Rb, 96 Rb, 97 S- 98t 

"Y, 100 Y, 101 Zr, 102 Nb, 103 Zr, 104 Nb, 105 Mo, 106 Nb, 

107 Mo, 108 Tc, 109 Mo, 110 Tc, m Ru, 112 Rh, 113 Ru, 114 Rh, 

115 Rh, 116 Rh, 117 Pd, 118 Rh, 119 Pd, 120 Ag, 121 Cd and 
122 / 



^Ag. 
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First of all, the fission path in our five-dimensional 
configuration space must be supplied, that is a depen- 
dence between all generalized coordinates. This tra- 
jectory starts from the ground-state of the system and 
reaches the exit point of the barrier. The ground-state 
corresponds to the minimal deformation energy in the 
first well. The adiabatic barrier in the multidimensional 
configuration space is determined by using the least ac- 
tion principle [23|]. Therefore, a trajectory in the multidi- 
mensional space is obtained by minimizing the functional 
(U} by using a numerical procedure as in Ref. |15j . The 
deformation energy landscape, minimized versus the pa- 
rameters Si (£=1,2) and ai/a,2, is plotted in Fig. [T]as a 
function of the necking coordinate C and the elongation 
parameter R. In the lower panel of Fig. Q]the least action 
trajectory is also plotted. In Fig. [3J the adiabatic fission 
barrier V, the effective mass M and the estimated mass 
number Ai are plotted along the fission path. According 
to Eq. ([1]) the inertia is computed within the cranking 
approximation |24j using the expression 

3=1 *=1 

where M qiqj are the elements of the inertia tensor com- 
puted for the generalized coordinates and qj . The sec- 
ond barrier top is located at i?=ll fm and corresponds to 
a mass ~ 120. For a constant charge density, this ratio 
of the mass asymmetry addresses a heavy fragment with 
A\ — Z\=81 and Z\=51, these values being close to magic 
numbers. Thus, the second saddle point corresponds to 
a partition that includes a double magic fragment. How- 
ever, as mentioned in Refs. (27| . microscopic approaches 
to fission f28-30| established that the second saddle point 
is asymmetrical with a value compatible with the ob- 
served mass ratio of the fragment distribution. There- 
fore, the dynamical saddle configuration obtained within 
our model is checked by minimizing statically the defor- 
mation energy around i?=ll fm. The eccentricity £\ and 
the neck parameter C are kept constant. The detailed 
region is displayed in Fig. [3] confirming that the saddle 
point is located at A^ «120 and R wll fm. 

Now we are interested in determining the fission barri- 
ers that address all the analyzed partitions. The asymp- 
totic deformations of the two fragments are taken from 
the literature [2||. Thus, the shapes of the initial nucleus 
up to the second barrier and those of the final fragments 
are known. In order to avoid a complicated determi- 
nation of the minimal action path for each partition, a 
linear variation from initial values of the generalized co- 
ordinates Ei (£=1,2) and a±/a2 is postulated starting from 
the saddle of the second barrier configuration up to the 
final ones, characterizing the fragments at the end of the 
fission process. 

By assuming a final mass asymmetry the values of the 
deformation energy V and the inertia M are plotted in 
Fig. Has a function of the elongation R and the light 
fragment mass A2- Some general features of the fission 



barrier can be extracted. The shapes of the external bar- 
riers change dramatically upon the mass-asymmetry. For 
A 2 <110 partitions, a double humped barrier ocurs while 
for more symmetric channels (A% > 110) a triple humped 
barrier is obtained. Moreover, the exit point from the fis- 
sion barrier is located at lower values of the elongation 
R for A2 < 110 than for more symmetric channels. The 
inertia tends to increase when the mass asymmetry in- 
creases. It is clear that the inertia favors the production 
of a symmetric partition, that is A 2 > 110. 

The spontaneous fission half-life is inversely propor- 
tional to the zero point vibration energy and to the prob- 
ability of penetration through the fission barrier. If the 
zero point vibration energy is considered the same for 
all partitions, then the yields are proportional with the 
penetrability in each channel. 

With these ingredients, the calculated penetrabilities 
at zero excitation energies through the barrier are com- 
pared to experimental yields in Fig. [5] We implicitly 
assume that the penetration of the inner barrier is the 
same for all partitions, so that differences in the barrier 
transmission between channels are induced only by the 
external barrier. 

We obtain a very good agreement between the the- 
oretical penetrability distribution and the experimental 
yields for A2 < 110. The maximum theoretical value is 
at ^2 = 107, while the maximum experimental yield is at 
the same value. A sudden drop of theoretical penetrabili- 
ties is theoretically obtained for channels with A2 > 110. 
As previously noticed, this behavior is connected with 
the ground state shapes of fragments that become oblate 
for these channels. Perhaps these oblate shapes are not 
the best configurations during the penetration of the bar- 
rier, and the final ground state oblate configurations are 
obtained only after the exit from the fission barrier. 

Concluding, we computed the cold fission path in the 
potential energy surface of 252 Cf by using the two center 
shell model, based on the idea of the cold rearrangements 
of nucleons during the cold fission process. We obtained a 
satisfactory agreement with experimental yields, by con- 
sidering variable mass and charge asymmetry beyond the 
first barrier of the potential surface. We can see that the 
mass asymmetry changed from a symmetric to the asym- 
metric configuration in the vecinity of the second barrier, 
due to the influence of the magic numbers Z=50 and 
N=82, i.e. 4 20 Ag+5f 2 Sb partition. It was shown that a 
good agreement with experimental data can be obtained 
only if the fission path proceeds through this saddle con- 
figuration. The final transmissions depend only on the 
external barrier and strongly depend upon the final frag- 
ment shapes. Due to the fact that prolate shapes are 
more favorable for the fission process, the maximal values 
of the yields are shifted towards channels characterized 
by lower values of A2. 

Thus, the cold fission process of 252 Cf can be called 
Sn-like radioactivity, similar with the Pb-like radioac- 
tivity, corresponding to various heavy cluster emission 
processes. We call all these processes shortly n (cluster) 
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decays. The peak in the final distribution corresponds This work was supported by the contracts IDEI-119 

to 107 Mo, due to the mass asymmetry degree of freedom, and PN 09 37 01 02 of the Romanian Ministry of Educa- 

allowing a lower barrier from Sn to this nucleus. tion and Research. 
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